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A PROBABILITY ANALYSIS OF THE METHOROLOGICAL 
FACTORS CONDUCIVE TO ATRCRAFT ICING 
IN THE UNITED STATES 


By William Lewis! and Norman В. Bergrun? 
SUMMARY 


Meteorological icing data obtained in flight in the United States 
are analyzed statistically and methods are developed for the determina- 
tion of: (1) the various simultaneous combinations of the three basic 
icing parameters (liquid-water content, drop diameter, and temperature) 
which would have equal probability of being exceeded in flight in any 
random icing encounter; and (2) the probability of exceeding any speci- 
fied group of values of liquid-water content associated simultaneously 
with temperature and drop-diameter values lying within specified ranges. 
The methods are particularly useful in the design of anti-icing equip- 
ment intended to operate through the United States, to define simulta- 
neous combinations of the meteorological variables which could be 
encountered, and to ascertain the effectiveness of the equipment in 
withstanding the natural icing conditions to which it may be subjected. 
In addition, a mathematical basis is provided for the future statistical 
analysis of meteorological icing data that might be obtained throughout 
the world. 


INTRODUCTION 


The program of research in aircraft ice prevention which has been 
conducted by the NACA during the past several years has been directed 
primarily toward the development of practical methods for the design of 
thermal ice-prevention equipment for various airplane components. Since 
a rational ice-prevention design requires a knowledge of the physical 
characteristics of icing conditions, an important phase of the research 
program has been an investigation of the meteorological conditions соп- 
dueive to icing. 


The severity of an encounter with icing conditions is determined 
principally by four factors; namely, the liquid-water content, the 
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diameter of the drops, the temperature, and the horizontal extent of the 
conditions. The meteorological investigation has, therefore, been con- 
cerned with obtaining measurements of these quantities in a wide variety 
of natural icing conditions in order to establish the range and rela- 
tive frequency of occurrence of various values, and combinations of 
values, of these quantities. 


During the period 1945 through 1948, a considerable amount of data 
was collected by the Ames and Lewis Laboratories. These results, and 
discussions of various phases of the investigation, have been reported 
in references 1 through |. All of these results, together with similar 
data from observations made by other organizations, were used as a basis 
for a listing of estimated maximum icing conditions recommended for con- 
sideration in the design of anti-icing equipment (reference 5). The 
probable maximum values listed in reference 5 were estimated on the 
basis of what was considered a reasonable extrapolation from a limited 
amount of data, and the concept of the probability of encountering 
various simultaneous combinations of values entcred only in a subjective 
and gualitative way. 


In this report, the available data obtained in icing flights in the 
United States are analyzed statistically and the results are presented 
graphically in two forms. In the first form, the results are plotted to 
show the various combinations of liquid-water content, drop diameter, 
and ambient-air temperature which have equal probability of being 
exceeded. In the second form, plots are presented to determine the 
probability of exceeding any specified value of liquid-water content 
under the condition that the value is associated simultaneously with 
values of temperature and drop diameter lying within specified intervals. 
Both plots are based on preselected values of horizontal extent; however, 
methods for estimating values applicable to other horizontal extents are 
also presented. 


Although the results of this report are directly applicable only to 
the United States, an estimation of the conditions prevailing in other 
localities of similar climatic conditions can be made based upon the 
results presented herein. As meteorological icing data are accumulated 
for regions outside the United States, the methods of this report should 
provide a framework for placing the data on a statistical basis that is 
not limited in scope to the United States. 


SYMBOLS 
The following symbols are used throughout this report: 


A,B,C designation applied to three events occurring simultaneously 
but not necessarily independently 
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rj 


М . 


Pa 


chord length, feet 
mean-effective drop diameter, microns 


the smallest drop diameter for which drops will impinge upon 
a thermal element 


interval of drop diameter, microns 
Naperian base of logarithms 
collection efficiency of an airfoil section, percent 


factor by which liquid-water-content values are multiplied 
to make them applicable to other distances of cloud hori- 
zontal extent, dimensionless 


weight rate of water intercepted, pounds per hour per foot 
of airfoil span 


an integer 


an icing condition considered to be defined by specifications 
of the three variables, liquid-water content, temperature 
depression below freezing, and drop diameter 


probability of overloading any given element of a thermal 
anti-icing system 


probability of the occurrence of event А 


probability of the simultaneous occurrence of three 
events, À, B, and C 


probability of the occurrence of event B, under the condi- 
tion that event A will occur 


probability of the occurrence of event C under the condi- 
tion that both events А and B will occur 


probability that any random icing encounter will be char- 
acterized by a combination of values of the variables 
(liquid-water content, temperature depression below freez- 
ing, and drop diameter), in which each variable must, гев- 
pectively, equal or exceed a specified value of that 
variable 


probability defined by the Gumbel probability distribution 
and expressed by equation (АЈ) 


R(W,T') 


Е! 


R' (W,D,T') 


S 


S(D,W,T') 


St 
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probability that, for a random icing encounter, the following 
three conditions will simultaneously be fulfilled: 
(1) liquid-water content exceeding a specified value; 
(2) temperature within a specified interval; and (3) drop 
diameter within a specified interval 


probability thet drop diameter will have 8 value lying within 
а specified range (AD) 


probability that temperature will have а value lying within 
а specified range (AT) 


probability that liquid-water content will be equal to or 
greater than Wj under the condition that the temperature 
lies within a specified range (AT) 


probability defined by the probability function Q(T') for 
a fixed value of Т! 


probability function, determined by the variation with tem- 
perature depression below freezing for the most severe 
icing conditions, of the percent of icing encounters in 
which the temperature depression below freezing is less 
than certain specified values 


probability that the maximum liquid-water content in an icing 
encounter is equal to or greater than W, under the condi- 
tion that the temperature depression below freezing is 
equal to or greater than T! 


probability function expressing В in terms of М and Т! 


probability that the liquid-water content associated with 
the maximum value of drop diameter in an icing encounter is 
equal to or greater than W, under the conditions that the 
temperature depression below freezing is equal to or 
greater than Т! and the maximum value of drop diameter 
is equal to or greater than D 


probability function expressing R' in terms of W, D, and T! 


probability that the drop diameter is equal to or greater 
than D, under the conditions that the temperature depres- 
sion below freezing is equal to or greater than T' and 
the liquid-water content is equal to or greater than W 


probability function expressing S in terms of D, W, and T: 


probability that the maximum value of drop diameter in an 
icing encounter is equal to or greater than D, under the 
condition that the temperature depression below freezing 
is equal to or greater than Т! 
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в: (Рут) 


gt '(D) 


probability function expressing 6! in terms of D and т! 


composite distribution function consisting of parts of each 
of the functions S(D,0,0') and S'(D,0) 


airfoil thickness, percent of chord length 

temperature existing in an icing condition, ОҒ 

temperature depression below freezing, equal to 32-T, F° 

temperature depression below freezing at which the heat out- 
put of the thermal element can bring the surface tempera- 
ture of the element just to freezing temperature, F° 

interval of temperature, OF 


interval of temperature depression below freezing, КО 


a constant which denotes the mode in the Gumbel probability 
equation 


airspeed, miles per hour 
liquid-water content, grams per cubic meter 


value of liquid-water content corresponding to a probability 
value of 0.05 defined by the Gumbel probability equation 


value of liquid-water content corresponding to a probability 
value of 0.63 defined by the Gumbel probability equation 


interval of liquid-water content, grams per cubic meter 
independent variable in the Gumbel probability equation 
а parameter which depicts the concentration of a frequency 


distribution about the mode in the Gumbel distribution 
equation 


Subscripts 


any selected icing condition 


a representative value of a variable 


о 
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ANALYSIS 


The data used in this analysis consist of 1038 sets of values of 
liquid-water content, drop diameter, and temperature obtained during 
252 encounters with icing conditions. Values of liquid-water content 
and drop diameter were measured by means of the rotating multicylinder 
method. All data available at the time of preparation of this report 
were used in the analysis, including (1) data obtained by the Ames and 
Lewis Laboratories and published in references 1, 2, 3, 4, and 6; 

(2) data contained in the monthly reports of the Air Materiel Command 
Aeronautical Ice Research Laboratory, the major portion of which is pre- 
верфей in reference 7; and (3) data obtained by United Air Lines (refer- 
ence 8) and American Airlines (reference 9) during icing flights with 
Douglas DC-6 type airplanes. 


For the purposes of this analysis, the basic unit of data is the 
icing encounter which consists of all measurements made during flight 
through a single area of continuous or intermittent icing conditions. 
This area was separated from other icing conditions? by relatively 
large areas in which no icing was observed. The problem considered 
herein is to determine the probability that an icing encounter chosen 
at random will include conditions of a given severity averaged over а. 
certain distance. Since most of the rotating-cylinder observations 
represent averages over distances of about 3 miles in cumulus clouds 
and 10 miles in layer clouds, these distances are regarded for purposes 
of data reduction as standard values of horizontal extent, applicable 
to data from the two principal cloud types. The probabilities obtained 
on this basis may be applied to other distances by using the data con- 
cerning the relation between horizontal extent and average liquid-water 
content obtained from continuous records of the rotating-disk icing-rate 
meter. (See reference 3.) It has been assumed that the greatest values 
of liquid-water content and drop diameter measured by the rotating су1- 
inders during an icing encounter represent maximum values averaged over 
the standard distances of 3 or 10 miles. This assumption is probably 
approximately true, since a particular effort was made to obtain 
rotating-cylinder data during the periods of most rapid ice formation. 


Since the variations of temperature during a single icing encounter 
are usually not very great, the most severe icing conditions measured 
during a particular icing encounter then will be either the greatest 
Value of liquid-water content and the corresponding value of drop diameter, 
or possibly, the greatest value of drop diameter and the corresponding 
value of liquid-water content. Both of these possibilities for & severe 





SThe term "icing condition" as used herein denotes а state of the atmos- 
phere defined by а set of values of temperature, liquid-water content, 
drop diameter, and pressure altitude in which the temperature is below 
freezing and the liquid-water content is greater than zero. 
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icing condition are considered in the probability analysis which, for 
the sake of convenience, is developed in three phases. 


The first phase of the analysis is the classification of the avail- 
able data according to cloud type and geographical location; the second 
phase is the determination of the particular combinations of liquid- 
water content, drop diameter, and ambient-air temperature, for a given 
horizontal cloud extent, that have the same probability of being simul- 
taneously equaled or exceeded in a random icing encounter; and the third 
phase is the determination of the probability of exceeding any specified 
group of values of liquid-water content associated simultaneously with 
temperature and drop-diameter values lying within specified ranges. 


Classification of Data According to Cloud Type 
and Geographical Location 


For classification purposes in this report, clouds will be divided 
into two classes, cumulus clouds and layer clouds; and the area of the 
United States will be divided into three regions, the Pacific coast 
region, the plateau region, and eastern United States with boundaries 
as indicated in figure 1. This system of classification divides the 
observations into six cases as indicated in table I, which includes 
information relating to the number of icing encounters and the number 
of individual measurements of liquid-water content and drop diameter for 
each group. The two cases for which the greatest amount of data are 
available, eastern layer clouds (case 5) and Pacific coast cumulus 
clouds (case 2), are treated in considerable detail, The other cases 
are treated in as much detail as the data permit, with the exception of 
case 6 (cumulus clouds in eastern United States) which is omitted because 
of insufficient data. 


Determination of Particular Combinations of Liquid-Water 
Content, Drop Diameter, and Temperature Depression 
Having the Same Probability of Being Simul- 
taneously Equaled or Exceeded in a 
Random Icing Encounter 


Exceedance probability.- The three fundamental icing variables used 
to define an arbitrary lcing condition are liquid-water content W, 
drop diameter D, and temperature T. For purposes of probability anal- 
ysis, however, it is more convenient to refer temperature to the freezing 
point and to use temperature depression below freezing T' as a variable 
to substitute for the temperature Т. Thus, any particular icing condi- 
tion can be represented by selected values of the variables Мі» Di, 
and Т.!. The probability that any random icing encounter of standard 
extent will provide simultaneous values of each of the icing variables 
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equal to or greater than each of the selected values (i.e., W2W,, D2D,, 
and T'2T,') will be termed the "exceedance probability" Pe correspond- 
ing to the icing condition Wi, Di, and Tife 


The exceedance probability represents the probability of the simul- 
taneous occurrence of three events which are not necessarily independent. 
In order to evaluate Ре, use is made of the theorem of compound proba- 
bility (reference 10) which states that the probability of the simulta- 
neous occurrence of three events, A, B, and C, is equal to the probability 
of the occurrence of A times the conditioned probability of the occur- 
rence of В, under the condition that A will occur, times the condi- 
tioned probability of the occurrence of С, under the conditions that 
both A and B will oceur. This theorem can be expressed symbolically 
as follows: 


Papo = РА Ра(А) Po(A,B) (1) 


If it is considered that the probability of exceeding certain values 
of Т', М, and D is analogous to the probability of the occurrence of 
events A, B, and C expressed in equation (1), then the exceedante 
probability Ре may be represented by the product of three probability 
factors Q, R, and S, according to the equation 


Ре = QRS (2) 


In applying equation (2) to the meteorological variables, the prob- 
ability Q may be considered as a function of T' alone, R a function 
of Wand T', and В a function of D, T', and М. The functional rela- 
tions between Q, R, and S, and D, T', and W are denoted, respectively, 
by the expressions 6(1!), R(W,T') and S(D,T',W). The meaning of each 
expression is as follows: 


Q(T') a relation that defines the probability Q in terms of Tt 
and expresses the probability for а given value of Т! 
of а random icing encounter having a temperature depres- 
sion below freezing equal to or greater than Т! 


R(W,T!) а relation that defines the probability R in terms of W 
when T! is specified and expresses the probability for 
a given value of W of а random icing encounter having 
а liguid-water content equal to or greater than W, under 
the condition that the temperature depression below freez- 
ing is equal to or greater than Т! 


S(D,W,T!) а relation that defines the probability 8, in terms of D 
when W and Т! are specified and expresses the probabil- 
ity for a given value of D of a random icing encounter 
having а value of drop diameter equal to or greater 
than D, under the condition that the liquid-water content 
is equal to or greater than W and the temperature 
depression below freezing is equal to or greater than Т! 


2H 
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The physical significance of the term exceedance probability is 
best illustrated by graphical means. Figure 2 presents a three- 
dimensional plot with the variables W, D, and T' as the axes. Any 
icing condition of standard duration defined by selected values of these 
three variables, such as Wp, Пі, and T,', appears as a single point 04 
оп the plot. Corresponding to the point 04, there 18 an ехсеедапсе 
probability Pe, which expresses the probability that a random icing 
encounter will provide simultaneous values of each of the icing vari- 
ables equal to or greater than each of the selected values, М, D1, 
and Ту’. According to equation (2), the value of Pe; would be given 
by the equation 


Pe, = 41815 


Thus, when Pe is regarded as а constant, Ре. equation (2) defines a 
surface passing through 01, specified by all the combinations of И, T', 
and D which have the same probability of being simultaneously equaled 
or exceeded. This surface is shown as the equiprobability surface in 
figure 2, and represents the locus of all icing conditions which have 
the same probability of being equaled or exceeded. The quantitative 
value of exceedance probability applying to individual points lying on 
the surface is expressed by equation (2). 


It should be noted that since Q, R, and S represent probabilities, 
they cannot be greater than one nor less than zero. Moreover, Q is 
equal to one when Т! equals zero and continuously decreases ав Т! 
increases; R is equal to one when W is zero; and for any constant 
value of Т', В is а continuously decreasing function of М. Finally, 
S is equal to one when D is zero and, for any particular pair of 
values of T' and W, S is a continuously decreasing function of D. 
The point at which the equiprobability surface intersects the T' axis 
is determined by the value of Т! for which Q = Ре, since R and 5 
are both equal to one. Similarly, the surface intersects the 
Wand D axes at values of W and D which make В and 5, respectively, 
equal to Pe. А constant-temperature contour passing through the 
point 0) оп the equiprobability surface would be defined by the соп- 
dition T! = T,'!, from which it follows that Q = 91, and hence 
RS = Ре/9.. These values are depicted in figure 2. 


The preceding theory concerning exceedance probability must be 
modified slightly to make it applicable to the particular problem under 
consideration, in which the probabilities are to apply to the most 
severe icing conditions in an icing encounter. In most cases, the most 
severe icing condition is represented by the maximum value of liquid- 
water content and the corresponding (simultaneously observed) value of 
drop diameter. For the portion of the equiprobability surface which 
represents low values of liquid=water concentration and large values of 
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drop diameter, on the other hand, the critical condition is represented 
by the maximum value of drop size and the corresponding value of liquid- 
water content. To take account of the first of these possible combina- 
tions, the function R(W,T') was determined from the greatest values of 
liquid-water content for each icing encounter and the function S(D,W,T') 
was obtained from the corresponding values of drop diameter; and to take 
account of the second possible combinations, another function S'(D,T') 
was determined from the greatest values of drop diameter and a corre- 
sponding function: R'(W,D,T') was obtained from the corresponding values 
of liquid-water content. The functions R(W,T') and S(D,T') were then 
used to determine the main portion of the equiprobability surfaces, 
while the functions R'(W,D,T') and S'(D,T') were used for the portion 
representing large values of D and small values of W. 


Construction of the eguiprobability surfaces.- The procedure used 
in the eonstruction of equiprobsbility surfaces from а particular set of 
data was essentially that of determining graphically the relations des- 
cribing the functions Q(T'), R(W,T'), S(D,W,f'), R'(W,D,T!), and S'(D,T!). 
This procedure is illustrated in appendix A. 


The mathematical functions used to represent the various distribu- 
tions were chosen mainly on the basis of the empirical criterion of 
goodness of fit, although theoretical considerations also were given 
congiderable weight, particularly in the selection of & function to 
represent the distribution of maximum values in those cases in which 
each item in the distribution is the maximum of a group of observed 
values. The function chosen for representing such distributions is one 
originated by Dr. Е. J. Gumbel (reference 11) to describe a distribution 
of extreme values. This function will be referred to hereinafter as 
"Gumbel's distribution." 


Because of an almost complete lack of observational data in flight 
concerning icing conditions at very low temperatures, reliance has been 
placed on the results of laboratory experiments on the formation of ice 
crystals to aid further in the construction of the equiprobability charts. 
The experiments relied upon (summarized in reference 12) indicate that a 
critical temperature exists in the neighborhood of -40° Е, below which 
ice crystals are formed in very great numbers in atmospheric air when- 
ever conditions of saturation with respect to liquid water occur. On 
the basis of this fact, it is inferred that clouds at temperatures 
below -400 Е are composed either entirely or almost entirely of ice 
crystals; hence the liquid-water content is assumed to be zero at all 
values of Т!” greater than 72° Е. 


The methods outlined in the preceding paragraphs have been used to 
prepare constant-temperature contours defining the equiprobability sur- 
faces for values of P, equal to 0.1, 0.01, and 0.001 for cases 1 
through 5 (table I). These curves are presented in figures 3 through T. 
Details of the probability calculations for cases 2 and 5 are given as 
examples in appendix A. 


A 
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The equiprobability surfaces defined by the contours presented in 
figures 3 through 7 are constructed only for three exceedance probabil- 
ities, but it may be desired to ascertain a point which lies between 
the surfaces defined for the three different exceedance probabilities. 
If such a point is desired, it can be identified in T'-D-W space by 
logarithmic interpolation between the equiprobability surfaces, since 
the probability functions decrease approximately exponentially with 
inereasing values of Т!, W, and D. 


The inclusion of horizontal extent as an additional variable.- 
Experience has shown that the more localized a cloud formation, the 
higher will be its liquid-water content. The results of the foregoing 
analysis (figs. 3 through 7), however, are directly applicable only to 


conditions not differing greatly in horizontal extent from the average 


distance over which the measurements were made, namely, 3 miles in 
cumulus clouds and 10 miles in layer clouds. The results may be made 
applicable, though, to other values of horizontal extent by making use 
of the data contained in reference 3 concerning the relation between 
maximum average liquid-water content and distance along the flight path. 


In order to modify values of liquid-water content obtained for one 
horizontal extent to be applicable to another, it has been found neces- 
sary to make three assumptions. These assumptions are that (1) the 
variables, drop diameter, and temperature are assumed not to be signif- 
icantly related to horizontal extent; (2) the variation of average 
liquid-water content with horizontal extent is assumed to be independent 
of drop size and temperature; and (3) the results from reference 3 are 
assumed to be representative of conditions as they would be encountered 
in normal flight operations. This last assumption is uncertain because 
of the limited number of flights analyzed in reference 3 (11 flights in 
layer-type clouds and 26 in cumulus clouds) and because the flight paths 
were frequently chosen with the aim of prolonging and maximizing the 
exposure to icing conditions. Оп the basis of these assumptions, equi- 
probability surfaces applicable to other values of horizontal extent may 
be constructed by multiplying all values of liquid-water content from 
the probability curves of figure 3 to 7 inclusive by the variable fac- 
tor F presented in figure 8. 


From figure 8, it can be seen that the horizontal extent selected 
as being applicable to a particular icing encounter can have a consid- 
erable effect on the factor Е used to modify liquid-water-content 
values. This effect is particularly pronounced in the case of layer 
clouds where the horizontal extent can be very large. For example, if 
an icing encounter in layer clouds is considered to extend over 100 miles 
instead of 10, the factor to use to modify liquid-water-content values - 
is about 0.4 instead of 1.0. Consequently, it is necessary to consider 
the appropriateness for design purposes of the standard values of cloud 
horizontal extent (3 miles for cumulus and 10 miles for layer clouds). 
For the "instantaneous" and "intermittent" classifications of reference 5 
(classes I and II ав defined therein), the horizontal extents 
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of 1/2 and 3 miles, respectively, as assigned in reference 5, are con- 
sidered to be reasonable values for design purposes. For the "contin- 
uous" classification (class III), a horizontal extent of 20 miles 
appears appropriate for design. Although it is realized that contin- 
uous icing might well extend beyond 20 miles, reduction of the liquid- 
water content by a factor less than 0.8 (an approximate value of F 
corresponding to 20 miles horizontal extent) does not appear justified 
for design purposes. It should be noted that these three design values 
of horizontal cloud extent, shown in figure 8 by broken lines, are sug- 
gested values only and do not preclude the use of the curves for the 
estimation of the factor F by selecting other values of horizontal 
extent. 


The Determination of the Probability of Exceeding Any Specified 
Group of Values of Liquid-Water Content Associated Simul- 
taneously With Temperature and Drop-Diameter Values 
Lying Within Specified Ranges 


The equiprobability surfaces discussed in the preceding section 
provide information concerning the simultaneous combinations of liquid- 
water content, drop diameter, and temperature having the same probability 
of occurrence. Another useful expedient for the designer, however, would 
be a means for determining the probability of encountering any icing ^ 
condition not included by the values of liquid-water content, drop diam- 
eter, and temperature specified as design criteria for a particular com- 
ponent of a thermal system. The basis for such a technique lies in 
determining the probability of equaling or exceeding certain values of 
liquid-water content specified for a number of icing conditions. 


Probability of encountering icing conditions which exceed certain 
specified values.- If, for a particular element of a thermal system and 
a chosen value of horizontal cloud extent, combination values of liquid- 
water content, drop diameter, and temperature are specified which the 
element can just tolerate, the values would define a surface in 
T'-D-W space. Such a surface can be viewed as representing the locus 
of liquid-water-content values to which the element is critical over 
specified ranges of temperatures and drop diameters. Individual values 
on such a surface would be obtained through knowledge of the area, rate, 
and distribution of water-drop impingement for the particular element. 


For purposes of explanation, figure 9 is presented to show the 
locus of critical values of liquid-water content for an assumed element 
of a thermal system. The locus of liquid-water-content values will be 
termed a "marginal" surface to indicate that values of М lying above 
the surface cannot be tolerated by the element of the thermal system 
under consideration. In general, the shape of the marginal surface in 
the W-D plane for the case of а wing element having a given heat out- я 
put: шау be assumed to have а shape resembling а hyperbola. А hyperbolic 
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shape is appropriate because for drop diameters less than the critical 
diameter* for the element under consideration, an infinitely large 
amount of liquid water is tolerable; and, as drop diameter progressively 
increases, decreasing amounts of liquid water are tolerable. Іп any 
W-T' plane, however, W will be assumed to vary from the values 

of W at T' =0 to W =0 at some value of T,', say Ter', where 
the heating capacity of the element can bring the surface temperature 

of the element just up to freezing temperature (a clear-air condition). 


For convenience, the base of the marginal surface (the T'-D plane 
at W = 0) is divided into n rectangles of dimension AD and AT' such 
that variations are small in values of liquid-water content included by 
the projection of each rectangle on the marginal surface. The value of 
liquid-water content considered to be representative for each rectangle 
is designated as Wi іп figure 9. The values of Т!, T, and D corre- 
sponding to М) shall be called Ty', Ty, and Dy, respectively. For 
every set of W4-D4-T4' values, there exists a probability function AP; 
which expresses the probability that, for a random icing encounter, the 
liquid-water content will exceed Wi and the drop diameter and tempera- 
ture depression will be within the ranges AD and AT', respectively. 


The n values of AP; corresponding to the n rectangles in the 
T'-D plane represent the probabilities of the occurrence of an event 
(the overloading of the anti-icing element) in n different, mutually 
exclusive ways. The probability that the event will occur is the sum 
of the probabilities that it will occur in any number of different ways 
provided the different possibilities are mutually exclusive. (See 
reference 10.) Hence, the probability of overloading the given anti- 
icing element is approximately equal to the sum of the partial probabil- 
ities for the n rectangles according to the equation 


n 
Р = > ДР: (3) 


In order to evaluate AP; for each of the п rectangles, it 
would appear that a separate determination of the relationship 
between W and AP; for each rectangle would be necessary. This pro- 
cedure is not practicable, however, because if the number of rectangles 
is chosen large enough to make variations in W. small for each rec- 
tangle, then the amount of data available to define AP; as a function 
of W is too small to yield statistically reliable results for the 
individual rectangles. Moreover, the presentation of the results would 
require separate curves or tables for each rectangle. It is possible 
to overcome these ditficulties by utilizing the fact that the 





“Тһе critical diameter may be defined as the smallest diameter for which 
drops will impinge on the element. For some airfoils, the critical 
size can be estimated by a technique outlined in NACA TN 2476 (refer- 
ence 13). и 
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distribution of drop diameter is approximately independent of tempera- 
ture and liquid-water content. Although there is a slight tendency 
toward higher frequencies of large values of drop diameter when the 
liquid-water content is low, this effect is too small to have an impor- 
tant influence in the determination of the total marginal probability. 
By making the assumption that the distribution of drop diameter is 
independent of temperature and liquid-water content, it is possible to 
determine values of AP; for different combinations of temperature 
depression below freezing, liquid-water content, drop diameter, and 
horizontal extent. 


When, for any given cloud horizontal extent, drop-diameter distri- 
bution is assumed to be independent of temperature and liquid water, the 
probability equation 


AP; = Pap Pap Pg, (AT) (1) 


can be written from analogy with equation (1). The interpretation of 
equation (4) is that the partial probability АР; applicable to а 
chosen point on the marginal surface is the product of three probabil- 
ities: the probability, Pap, that the drop diameter will be within the 
interval AD; the probability, Pam, that the temperature will lie within 
the interval АТ; and the conditioned probability, Py, (AT), that the 
liquid-water content will be equal to or greater than Мі under the 
condition that the temperature lies within the interval AT. 


To evaluate equation (№), values of Pap were obtained from the 
frequency distribution of all observations of drop diameter; values 
of Par were obtained from the frequency distribution of the tempera- 
ture of the most severe? icing observations per icing encounter; and 
values of Py, (AT) were obtained from the frequency distributions of 


maximum liquid-water content per encounter for the various temperature 
intervals. Because of the assumption that drop size is independent of 
temperature and liquid-water content, it was not possible to provide 
exactly for the case in which the maximum observed drop diameter deter- 
mines the most severe icing condition in an icing encounter. An approx- 
imate allowance was made for this factor, however, by substituting the 
frequency distribution of the greatest value of drop diameter per 
encounter over & limited range of large drop diameters. 


Equation (4) was evaluated for different horizontal extents, and to 
do this the frequency distribution for maximum liquid-water content per 
encounter was altered according to а relation from reference 3 between 
maximum average liquid-water content per encounter and the distance 
along the flight path. 


ee 
“For the purposes of this report, the most severe icing condition is con- 
sidered as the one having the largest liquid-water content. 
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The partial probability charts.- In evaluating equation (4) for 
any one value of cloud horizontal extent, numerous values of AP; are 
possible from the meteorological data for different combinations 
of Pap, Par, and Py, (AT). The values of ДР4, however, сап be conven- 
lently summarized in the form of partial probability charts. Such 
charts are presented in figure 10 for cases 1 through 5 of table I. 


As is evident from figure 10, seven intervals of drop diameter and 
six intervals of temperature were used in the construction of the 
charts. These intervals were chosen so as to avoid large variations 
of W encompassed by any particular AD-AT' rectangle projected on 
the marginal surface. Ти all conceivable combinations, the temperature 
and drop-diameter intervals provide 42 different values of AP; which 
must be summed (in accordance with equation (3)) to obtain the probabil- 
ity, P, of overloading the particular element. 


To obtain a specific value of АР) from a chart of figure 10, the 
following procedure is used: Choose a value of liquid-water content Wy 
and follow a vertical line until it intersects with the desired 
horizontal-extent line; from this intersection, follow horizontally until 
an intersection with the appropriate temperature-interval curve is 
obtained; from the point of this intersection, follow a vertical line 
downward to obtain an intersection with a chosen drop-diameter curve; 

а horizontal line through the latter intersectian yields the desired 
value of АР: on the ordinate scale. For ease in observing this pro- 
cedure, each chart in figure 10 has a system of arrows leading from 
some arbitrary value of Wi through the horizontal-extent, temperature, 
and drop-diameter curves to a corresponding value of АР. 


It should be noted that partial probabilities presented in figure 10 
can be determined, for values of horizontal extent not specifically 
represented in the charts, by interpolation between the horizontal-extent 
сигуев. It is not permissible, however, to interpolate between the 
curves representing temperature and drop diameter since the partial prob- 
abilities are dependent wpon the size of the intervals for which the 
temperature and diameter curves are constructed. Тһе approximate median 
values of temperature and drop diameter, which correspond to the intervals 
of Ty and Пі used in the charts of figure 10, are listed in table ТТ 
for reference. Median values are used because the median has the advan- 
tage that its value is not greatly influenced by the chance variations 
of extreme items. A detailed description of the techniques used in 
constructing the partial probability charts, with application to selected 
cases, is presented in appendix B. 
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DISCUSSION 


Reliability of the Equiprobability and the 
Partial Probability Charts 


The reliability of the probability analysis reported herein is 
affected chiefly by four factors; namely, (1) errors of measurement, 
(2) amount of data available, (3) representativeness of the data, and 
(4) limitations with respect to climate and altitude. The possible 
effect of these factors is discussed in the following paragraphs. 


Errors of measurement.- As noted in reference 14, the measurement 
of liquid-water content and drop diameter by the multicylinder method 
ін subject to error. The following factors may be listed as possible 
sources of error: (1) runoff, (2) bounce-off, (3) evaporation, 
(4) variations in drop concentration in the neighborhood of the cylinders 
due to disturbing influences of the airplane, (5) collection of snow 
flakes as well as liquid drops in mixed clouds, (6) inaccuracies due to 
the theoretical drop-size distributions differing from those assumed in 
the calculations, (7) failure to fulfill, with the apparatus used, the 
theoretical assumption of two-dimensional flow, (8) irregularities in 
the surfaces of the ice-covered cylinders, (9) errors in weighing, 
(10) irregular breaking of the ice when separating cylinders, (11) errors 
in measuring the final diameter of the ice-covered cylinders or in cal- 
culating the average diameter from an assumed value of ice density, 
(12) errors in the measurement of airspeed and exposure time, and 
(13) errors in matching the experimental data with the theoretical 
collection-efficiency curves. 


A detailed discussion of the possible effects of each of these 
factors on the accuracy of the probability analysis presented herein is 
beyond the scope of this report. The problem of runoff, though, is 
worthy of discussion since this factor imposes an upper limit to the 
range of liquid-water-content values that can be measured. The results 
of heat-transfer calculations reported in reference 15 indicate that the 
maximum value of liquid-water content that can be measured with rotating 
cylinders as they are normally used varies approximately linearly with 
temperature depression below freezing, reaching about 2 grams per cubic 
meter at 5° Е for an airspeed of about 200 miles per hour, An examina- 
tion of the data used in the present analysis indicates that not more than 
about 5 percent of the icing encounters in any of the cases analyzed соп- 
tained observed values of liquid-water content which were close enough to 
the theoretical maximums to be regarded as likely to be influenced by 
runoff. Ти order to estimate the effect of a small percentage of errors 
of this type on the results of the probability analysis, it is necessary 
to consider the manner in which the liquid-water-content distribution 
funetions were determined from the experimental data. Ти this process, 
the cumulative frequency distribution of maximum liquid-water content 
per icing encounter was plotted on Gumbel's distribution paper, and the 
straight line of best fit was drawn to represent the entire distribution; 
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m thus errors in the upper 5 percent of the distribution had only а minor 
effect in the final selection of the line used to represent the distri- 
bution function. Moreover, it was noted that in nearly all cases the 
entire distribution followed the straight line with minor and apparently 
random variations and without any consistent or marked tendency for large 
negative departures at the upper end of the distribution. It may be 
inferred, therefore, that the limitation on the measurement of liquid- 
water content discussed in reference 15 does not have an important influ- 
ence on the results of this analysis. 


In view of the number and complexity of the possible sources of 
error, it is not possible at this time to make a reliable estimate of the 
total accuracy of the measurement of liquid~water content. In а statis~ 
tical analysis, however, reasonably reliable results may be obtained as 
long as the errors are small compared to the real variations of the 
quantity measured and are distributed approximately at random. It is 
believed that these conditions are approximately fulfilled for the values 
of liquid-water content measured by the multicylinder technique. 


Multicylinder measurements of drop diameter, on the other hand, while 
quite accurate and reliable for small drops, become increasingly inaccurate 
as the drop size increases. Moreover, the errors are not normally dis- 
tributed, since large positive errors are more probable than large nega- 
tive errors, especially at large values of drop diameter (see reference 1h). 
The result of these errors is an increase in the dispersion of the fre- 
quency distribution and an exaggeration of the probability of occurrence of 
к large values of drop diameter. This effect is probably not important for 

observations in the eastern United States because in this case the drop 
size rarely exceeds the range of reliable measurement. In the Pacific 
coast area, on the other hand, where large drops occur with much greater 
frequency, it is quite likely that the extreme values of drop diameter are 
influenced considerably by this factor. 


Amount of available data.- Some concept of the accuracy of the ргоба- 
bility analyses for the various cases may be obtained by reference to the 
number of icing encounters for each geographical area and each cloud-type 
classification. These data are presented in table I. It is noted that 
the greatest amount of data for any one case is 110 encounters in layer 
clouds in the eastern United States. A sample of this size is sufficient 
to establish the probabilities with considerable confidence at the 0.01 
probability level, and the extrapolation to a probability of 0.001 is not 
likely to introduce serious errors. In the remaining cases, the sample 
size ranges from 25 to 44 encounters. In these cases, the extrapolation 
to a probability of 0.001 may introduce considerable uncertainty. 


Š Representativeness of the data.- Another factor, which probably has 
as great an influence on the reliability and applicability of the analy- 
sis as either errors of measurement or sample size, is the systematic bias 

~ introduced by the fact that all the test data were obtained during flights 
in which icing conditions were deliberately sought instead of being aeci- 
dentally encountered. It is difficult to estimate the effect of this 


18 МАСА TN 2738 


factor since measurements of liquid-water content and cloud-drop size 
encountered during ordinary transport operations are not yet available. 

It would appear that the principal difference between transport experi- 
ence and icing research flights would be the greater frequency of encoun- 
ters with icing conditions in the latter. This difference alone, however, 
would not affect the validity of the analysis since the icing encounter is 
taken as the unit of experience upon which the analysis is based. The 
question to be considered is whether the icing encounters which occurred 
during the research flights represent an unbiased sample of the encoun- 
ters which wouid occur during normal flight operations. 


The differences in severity and extent between icing conditions 
experienced during research flights and the conditions which would have 
been encountered during an ordinary flight in the same general area at 
the same time are of two types. The first type occurred in cases in 
which no icing would have been encountered on a normal flight through 
the area but, due to the existence of clouds in thin layers or localized 
areas, it was possible to experience icing with the test airplane. 
Encounters of this type would be expected to be less seyere and exten- 
sive, on the average, than typical encounters during normal operations. 
Differences of the second type occurred in cases in which icing condi- 
tions would һауе been encountered on а normal flight but, due to а con- 
scious effort on the part of the pilot, the conditions experienced by 
the research airplane were more severe and extensive. It is evident, 
therefore, that the special procedures followed during research flights 
have resulted in the inclusion in the test results of data from icing 
encounters which were both more and less severe than those which would 
have been experienced under normal operating conditions. It is not 
possible to state with certainty what is the net effect of these two 
types of differences, but it is believed that differences of the second 
type predominate, with the result that the data are somewhat, though not 
greatly, biased by the inclusion of an abnormally high frequency of 
severe lcing conditions. 


Effect of climate and altitude.- Another factor having an important 
bearing on the applicability of the results is the effect of climate and 
flight altitude on the temperature of icing conditions. It is likely 
that the distribution of temperatures of icing conditions encountered by 
the test airplanes closely resembles that which would be found in normal 
operations of unpressurized airplanes during winter and spring in the 
areas studied. Different temperature distributions would be expected, 
however, With high-altitude airplanes or with conventional airplanes in 
other climates and seasons. For example, if an airplane is expected to 
cruise at altitudes of from 18,000 to 28,000 feet, where the temperature 
is normally in the range between O and -40° Е, a much higher frequency 
of icing at low temperatures is to be expected and the results of an 
analysis of icing conditions occurring mostly at higher temperatures 
would be applicable only to a limited extent. 
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Comparison of Results of the Equiprobability 
Analysis With Results of NACA TN 1855 


Certain values of liquid-water content, drop diameter, and temper- 
ature were presented in NACA TN 1855 (reference 5) for consideration in 
the design of anti-icing equipment. The combinations of meteorological 
variables recommended, however, were lacking in one respect, namely, the 
probability of occurrence of these particular combinations. Fortunately, 
with the aid of the equiprobability curves, some indication of these 
probabilities can be obtained when corresponding meteorological condi- 
tions are compared. 


Four classes of icing conditions are presented in NACA TN 1855: 
I, instantaneous; II, intermittent; III, continuous; and IV, freezing 
rain. All classes, except IV, are subdivided into two types of icing 
conditions, maximum and normal, and are confined to icing conditions 
associated with definite cloud formations of characteristic horizontal 
extents. Hence, classes I, II, and III are ideal for comparison with 
the results which can be obtained from the equiprobability charts 
(figs. 3 to 7) and the relations for horizontal extent (fig. 8). 


Two groups of icing conditions are chosen as a basis for comparison 
between the analysis of this report and NACA TN 1855. ‘The first group 
consists of icing conditions selected from the probability analysis for 
the case of layer-type clouds with an exceedance probability 
of Pe = 0.001. The second group consists of icing conditions for the 
case of Pacific coast cumulus, also with an exceedance probability 
of 0.001. These two groups are compared, respectively, with classes III-M 
and ТТ-М from NACA TN 1855. Тһе comparisons are made by determining 
values of liquid-water content from figure 4 and figures 3, 5, and 7 for 
the various combinations of temperature and drop diameter given in 
table I of NACA TN 1855 for class II-M and III-M conditions, respectively. 
The three values of liquid-water content obtained from figures 3, 5, 
and 7 corresponding to class ТТТ-М were averaged by weighting the values 
according to the horizontal extent of the region to which they apply. 

The three regions (fig. 1) have area ratios of about 1/8, 2/8, and 5/8 
for the Pacific coast, plateau, and eastern regions, respectively; and 
these area ratios were used to weight the values of liquid-water content 
from the equiprobability charts so that the resulting average values 
would be .оп a comparable basis with the values expressed in NACA TN 1855. 
The weighted average values of liquid-water content were corrected by 
figure 8 to make them applicable to a horizontal extent of 20 miles, a 
figure reasonably applicable to the layer-cloud data of NACA TN 1855. 
For class ТТ-М conditions, no correction for horizontal extent was 
required because the values in NACA TN 1855 and figure 4 both apply to 

3 miles horizontal-cloud extent. 


The results of the comparisons between NACA TN 1855 and the proba- 
bility analysis are shown in tables III and IV for classes II-M and III-M, 
respectively. Іп general, the liquid-water-content values agree quite 
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closely, except at temperatures of -MO F and below. Іп this low tempera- 
ture range, the values of liquid-water content derived from the equiprob- 
ability surfaces are considerably lower than values given in МАСА TN 1855. 
The reason is that the probability analysis entails some extrapolation in 
this temperature range, and so may not provide exactly the correct value 

ОҒ liquid-water content. Оп the other hand, the values listed in 

NACA TN 1855 are not restricted to actual measurements, and therefore may 
not be of precise magnitude. For design purposes, however, the values 
given by the probability analysis should be of proper magnitude. Іп this 
regard, it should be noted that the majority of the data utilized in the 
probability analysis was taken at comparatively low altitude (13,000 feet), 
whereas the temperature range between ® F and -40° Е represents consider- 
ably higher altitudes (18,000 to 28,000 feet). This difference in altitude 
could have some bearing on the accuracy of the probability data, particu- 
larly at very high altitudes. 


By a procedure similar to that used in the comparisons for the two 
example cases, the correspondence between classes of icing conditions 
represented by the equiprobability charts and other classes of conditions 
presented in NACA TN 1855 could be ascertained approximately. The согге- 
вропдепсе for these classes of conditions, and also for the two example 
cases, are shown in table У. 


An inspection of table V reveals a consistency between the, classes 
of design values recommended in NACA TN 1855 and specific values of 
exceedance probabilities. Only the instantaneous maximum icing condi- 
tion (Ре = 0.0001) appears to be incongruous with the maximum conditions 
in classes ІІ-М and ІІІ-М (Ре = 0.001). This apparent incongruity is 
directly attributable to the fact that the instantaneous maximum condi- 
tion presented in NACA TN 1855 was calculated for tall tropical cumulus 
clouds. Such extremely severe icing conditions would, of course, be 
exceeded infrequently as is borne out by the probability analysis. 
Another point of interest in regard to table У is the general order of 
magnitude of the exceedance probabilities which are found to apply to 
the various classes of icing conditions specified in NACA TN 1855. For 
example, any one of the icing conditions listed under class II-M (inter- 
mittent maximum) or class ІІІ-М (continuous maximum) will be exceeded once 
in only about a thousand encounters with icing conditions. It is believed 
that the icing conditions specified by classes II-M and 111-М do not 
impose too severe design requirements and, therefore, if the ice- 
prevention system will cope with these conditions, the results of the 
probability analysis indicate that satisfactory ice protection would be 
provided for the vast majority of icing encounters. 


USE OF THE EQUIPROBABILITY AND PARTIAL PROBABILITY 
CHARTS AS APPLIED TO THE DESIGN OF THERMAL ТСЕ- 
PREVENTION EQUIPMENT 


Heretofore, no information has been available regarding the prob- 
ability of encountering icing conditions which exceed certain specified 
values, so that some discussion of probability as applied to the design 
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of thermal anti-icing equipment seems worthwhile. Іп general, two design 
problems exist: one problem is to establish consistent simultaneous com- 
binations of liquid-water content, drop diameter, and temperature which 
can occur in nature, and the other problem is to determine how effectively 
a given thermal element of an anti-icing system will cope with the entire 
array of icing conditions existing in nature. It is the purpose of the 
equiprobability charts to aid in the solution of the first problem, while 
the partial probability charts are intended to aid in the solution of the 
second probiem. 


The Equiprobability Charts 


The equiprobability charts essentially are a means to determine 
quickly consistent sets of liquid-water content, temperature, and drop- 
diameter values which have the same probability of being exceeded during 
icing encounters. Sets of Buch values obtained from the equiprobability 
charts can be used as criteria for the design of thermal anti-icing 
equipment. No need to use these charts exists if it is found that 
table V will define a class of conditions from МАСА TN 1855 suitable for 
the needs of the particular design problem. However, should it be 
desired to design, вву, & jet-engine inlet duct for ал ехсеедалсе prob- 
ability of 0.001 (only one chance in а thousand of exceeding а specified 
amount of liquid water), table V would not be of help because correspond- 
ing icing conditions from NACA TN 1855 are not defined for this partic- 
ular case. Hence, reliance would have to be placed on the equiprobapility 
charts in order to establish consistent sets of values for a chosen 
value of Pa. 


General procedure for using the equiprobability charts,- The follow- 
ing procedure for the use of the equiprobability charts is suggested: 


1. Determine for what class of icing condition (I, II, or III of 
WACA TN 1855) a design is to be made by considering the nature of the 
component to be protected. (See table V.) | 


2. Establish the severity of icing condition within the class by 
selecting the exceedance probability for which the component is to be 
designed. (See table V.) 


3. With the aid of table V and steps 1 and 2, determine what equi- 
probability chart to use in defining values of W consistent with dif- 
ferent sets of Т and D. The selection of charts can be made by noting 
the type of cloud formation applicable to the class (i.e., I, II, or III) 
and severity of icing condition established in steps 1 and 2. For cumulus 
clouds, it is suggested that the charts for the Pacific coast area be 
used for design purposes rather than the charts for the plateau area. 

The reason is that, until more data are available for the eastern 
United States, Pacific coast cumulus-cloud data should provide liquid- 
water-content values that are conservative for the remainder of the 
United States. For layer clouds, however, it is suggested that the 
charts for all three geographical areas be used for design purposes by 
using a weighted average of liquid-water-content values from the three 
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charts. The justification for using a weighted average is that layer 
clouds can be regarded as being structurally continuous over exceedingly 
large geographical areas; consequently, only average values of the 
parameters defining an icing condition would seem to have much signif- 
icance in the design of an anti-icing component for layer clouds exist- 
ing throughout the United States. 


h. From the appropriate equiprobability chart, choose values 
of М corresponding to а variety of sets of T and D. Ти the case of 
layer clouds, where more than one chart may be employed to determine 
values of W, a suitable weighted average for the entire United States 
can be obtained by weighting the liguid-water-content values for the 
Pacific coast, plateau, and east coast regions, in the ratios of 1/8, 2/8 
and 5/8, respectively. These weight ratios, 1/8, 2/8, and 5/8, are pro- 
portional to the three geographical areas and may be considered as 
approximately accounting for the length of time that an airplane flying 
transcontinentally would be in each area. 


5. Multiply the values of М determined in step № by a factor, 
chosen from figure 8, for including the effect of horizontal extent of 
cloud formation. The factor chosen should lie within the range of 
horizontal extent indicated in figure 8 Рог the particular class icing 
condition established in step 1; it should correspond also to the curve 
for the exceedance probability decided upon in step 2. Recommended 
values of the factor correspond to values of horizontal extent identi- 
fied in figure Ө ав design values. These design values are for cloud 
horizontal extents of 1/2, 3, and 20 miles and correspond to icing con- 
dition classes I, II, and III, respectively. 


Example showing the use of the equiprobability charts.- For this 
example, a design of a jet-engine inlet-guide vane will be assumed in 
which it is desirable to protect against those icing conditions having 
a probability of being exceeded once in a thousand icing encounters 
(Ре = 0.001). The steps in determining sets of values of W, T, and D 
are? 


1. From table V, jet-engine intakes are recommended to be designed 
for class I icing conditions. 


2. The exteedance probability for which the design is to be made 
is Ре = 0.001. 


3. In table У, Pacific coast cumulus clouds are considered typical 
for class I icing conditions. Hence, equiprobability charts for 
Pacific coast cumulus (fig 4(c)) are chosen for obtaining values of 
liquid-water content corresponding to various combinations of temperature 
and drop diameter. 


h. Values of М selected from figure (c) for different arbitrar- 
ily chosen sets of T and D are presented in the following table: 
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Instantaneous Conditions Corresponding to 
[Horizontal extent not considered] 





5. From figure 8, it is determined that the factor by which to 
multiply the liquid-water-content values presented in step | is 1.3. 
This particular value corresponds to the design value of cloud horizontal 
extent suggested in figure 8 for class I (instantaneous) icing conditions 
with Ре taken as 0.001. A tabulation of the final liquid-water-content 
values suitable for design purposes is presented in the following table: 


Instantaneous Conditions Corresponding to 
Pe = 0.001 
[Horizontal extent considered] 
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In the above table, care should be taken to note that the exceedance 
probability, Ре = 0.001, applies to each icing condition listed and not 
to the class in general. 


The Partial Probability Charts 


The partial probability charts essentially provide a rapid means of 
determining the approximate probability of exceeding any specified group 
of values of liquid-water content associated simultaneously with tempera- 
ture and drop-diameter values lying within specified ranges. Thus, if 
an anti-icing component (e.g., unit span of а jet-inlet guide vane) is 
designed by the method of reference 16 to withstand a constant value? 
of weight rate of water-drop impingement Mt, as computed from а partic- 
ular set of values of liquid-water content, temperature, and drop dian- 
eter, and if other combinations of these same variables are found which 
correspond to the same weight rate of drop impingement, then the prob- 
ability of encountering conditions under which the component will not 
perform satisfactorily can be determined approximately. The calculation 
of this probability can be performed with the aid of the partial prob- 
ability charts (fig. 10). 


General procedure for using the partial probability charts.- To use 
the partial probability charts, the following procedure is suggested: 


l. Establish values of liquid-water content, Мі; for various com- 
binations of drop diameter and temperature, which are regarded as crit- 
ical (a marginal surface) for a given component of а thermal system. 
Values of Мі сап be obtained by computing values of liquid-water content 
required to provide a constant weight rate of water impingement on the 
component. The equation " 

t 


ИН) = ———— 
1 0.33E,Vte 


(5) 





Тре choice of a constant value of weight rate of water-drop impingement 
rather than, for example, the use of the amount of heat to raise the 
surface temperature to 32° F, provides, according to reference 16, a 
reasonable basis for estimating the probability of encountering an 
icing condition too severe for the ice-protection equipment. The 
amount of heat supplied for the particular weight rate of water impinge- 
ment selected depends upon the type of protection required for the 
component, that is, whether, for example, all the water impinging 
should be evaporated instead of having the surface temperature merely 
brought to 329 Е. Since for various water-drop sizes and temperatures 
the amount of heat required for protection is not unique for a given 
weight rate of water-drop impingement, the marginal surface will not 
represent the same degree of protection for all combinations of drop 
size and temperature. However, it will very nearly do so. As a result, 
the probability value computed from use of the marginal surface is not 
exactly correct. 
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derived in reference 16, can be used to calculate the values of Wi as 
а function of drop diameter and temperature, since E, is a function 
of these latter two variables. The particular value of M, used to 
calculate values of Wi is obtained by inserting into equation (5) 
initial values of Wi, Ет, Vy t, and с corresponding to а selected 
meteorological design condition and airfoil section. The appropriate 
values of E, 50 employ in equation (5) for the chosen values 

of V, t, and с must be determined by using some type of water-drop- 
trajectory calculation, such as presented in reference 13. When using 
equation (5) to obtain values of Wi, the computations should be made 
for all the possible combinations of T4 and D, listed in table II. 
In total, there are 42 such combinations. 


2. Determine the probability, ЛАР), of exceeding each of the values 
of liquid-water content, Wi, determined in step 1. To determine а 
specific value of AP, from & particular partial probability chart 
presented in figure 10, the procedure is as follows: Follow & vertical 
line representing the liquid-water content until the line intersects & 
horizontal-extent line. From the horizontal-extent line, follow hori- 
zontally until an intersection with а temperature-interval contour is 
obtained. From this intersection, follow a vertical line downward to 
obtain an intersection with а diameter-interval line. A horizontal line 
through this latter intersection yields the value of ДР; оп the ordi- 
nate scale. When it happens that & horizontal line drawn from the 
horizontal-extent curye will not intersect the desired temperature- 
interval line, all that can be said about the value of AP, for a given 
diameter interval is that it is smaller than the value the horizontal 
line would indicate when projected to the extreme left side of the раг- 
tial probability chart. Such values are usually small enough to be 
neglected. 


3. Establish the probability, P, of exceeding 211 the values of 
liquid-water content defined in step 1. This probability may be obtained 
by summing all the individual values of ЛАР), determined in step 2, іп 
accordance with equation (3). The number of AP; values which must be 
sumed is 42, since this number represents all combinations of tempera- 
ture and drop diameter possible by using the charts of figure 10. 


Example showing the use of the partial probability charts.- The 
example of the jet-engine inlet-guide vane will be continued to demon- 
strate the use of the partial probability charts for a vane section 
assumed to be heated to 100-percent chord. The starting point of this 
example will be taken with established sets of values of liquid-water 
content, drop diameter, and temperature that define а marginal surface 
for a unit span of the vane. These sets of values, which are presented 
in table VI, will be presumed to have been obtained by calculating 
values of Wi from equation (5), using a particular design icing condi- 
tion from the equiprobability charts. The particular icing condition 
presumed to have been selected for the computations is one taken from 
the example use of the equiprobability charts and defined by the following 
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values: W = 2.48 grams per cubic meter; T = 09 F; D = 15 microns; and 
cloud horizontal extent, 0.5 mile. With these available data, the steps 
employed to use the partial probability charts are: 


1. Values of AP; are selected from the marginal probability 
charts for Pacific coast cumulus (fig. 10(b)) to correspond with the 
conditions chosen from the equiprobability charts for a horizontal 
extent of 0.5 mile. The results are tabulated in table VI for the 42 
different possible combinations of Ті and Di. 


2. The probability of encountering &n icing condition lying out- 
side the surface defined in W-T'-D space by the W4-T4-D4 values 
listed in table VI is obtained by adding the partial probabilities found 
for the various conditions presented in the table. In the summation, 
values entered in table VI ав being less than a certain specified amount 
are neglected because the values are exceedingly small and not exactly 
determined in magnitude. These values arise from the fact that the 
charts are not constructed to encompass quite small values of АР}. For 
the example presented in table VI, the value of the probability, P, was 
found to be 0.0657. The value of 0.0657 can be interpreted to mean that 
in approximately 93 out of 100 icing encounters the values of liquid- 
water content, for all combinations of drop diameter, will lie below the 
marginal surface and hence be encompassed by the design. 


CONCLUDING REMARKS 


The equiprobability surfaces and marginal probability calculation 
charts presented herein provide, for the United States, a representation 
of the data now available on the meteorological factors responsible for 
aircraft icing, expressed in terms of the probability of exceeding 
various combinations of values of these factors. The usefulness of the 
results is limited somewhat by the amount of data available; errors of 
measurement, especially at large values of drop diameter; nonrepresent- 
ativeness of the flight procedures during research flights; and limited 
altitude and geographical extent of the research flight program. 


In spite of these limitations, the probability analysis presented 
does provide an indication, heretofore unavailable, of the combinations 
of icing conditions having equal probability of being exceeded in the 
United States, and also the probability of exceeding a special set of 
icing conditions. Also, a procedure has been established for the в%а- 
tistical analysis of future icing meteorological data obtained on a 
world-wide basis. 


Ames Aeronautical Laboratory 
National Advisory Committee for Aeronautics 
Moffett Field, California, April 18, 1952 
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APPENDIX A 
DETATLS OF DETERMINING THE EQUIPROBABILITY 


SURFACE CONTOURS 


The distribution functions used in the equiprobability study are 
the normal and the Gumbel probability distributions. Graphical methods 
were used to evaluate the distributions by using normal-probability 
graph paper for the normal distributions and a specially constructed 
` paper for the Gumbel distributions. These special types of graph paper 
have the property that the cumulative distribution turves are repre- 
sented by straight lines. 


The Gumbel distribution, taken from reference 11, is given by the 
equation 


Pa = 1-е-е"АХ-Ш) (A1) 
In equation (Al), the constant u is the mode and о is a quantity 
which measures the concentration of the frequency distribution about the 
mode. Since on the specially constructed Gumbel distribution graph paper 
equation (Al) will appear as a straight line, it is desirable to define 
the line by specifying two points on the curve. For the purposes of 
this analysis, the points are chosen at which the probability Pq has 
the values of 0.63 and 0.05. The value of Ра = 0.63 may be found by 
placing х =u and the value of Ра = 0.05 may be found by 

placing a(x-u) = 3.0. 


When a straight-line fit could not be obtained from the data by . 
using Gumbel distribution paper, an arbitrary smooth curve on normal 
probability paper was used to represent the distribution function. This 
procedure conserved the essential form of the observed distribution, 
smoothed out the irregularities, and provided a reasonable basis for 
limited extrapolation. 


Cases 2 and 5 (table I) are used to demonstrate how the normal and 
Gumbel probability distributions were applied to the data since these 
cases were the ones in which the largest number of encounters with icing 
were recorded. As mentioned in the text, the general procedure is to 
evaluate the functions Q(T'), R(W,T'), S(D,W,T'), R'(W,D,T') 
and S'(D,T'). 


Case 2, Cumulus Clouds in the Pacific Coast Area 
The distribution functions.- The distribution of the temperature 


depression below freezing of the most severe icing condition in each 
icing encounter for this case is shown in figure ll plotted on Gumbel's 
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distribution paper. Іп this figure the broken line represents the 
observed values and the straight line is the estimated line of best fit. 
This line defines the function Q(T!). 


The distribution of maximum liquid-water content per encounter was 
plotted for all values of Т' equal to or greater than each of the 
following values: 0, 8, 12, 17, 22, 27, 28, 32, and 37. The plots 
for Т!20, and T'>27 are presented in figure 12. The values of the 
distribution parameters, Ио вз and Wo „og, obtained from these distribu- 
tions and plotted as functions of T', are shown in figure 13. The 
smoothed curves in figure 13 define the liquid-water-content distribu- 
tion function, R(W,T'). These curves were drawn to the spontaneous 
freezing temperature, Т! = 72 FO. 


The distribution, S(D,W,T') of drop size to be used with the dis- 
tribution of maximum liquid-water content per encounter is to be deter- 
mined for all observations with Wand Т! greater than certain values. 
Small values of W need not be considered in this case since conditions 
of small М and large D are to be expressed in terms of 5! and R'. 
The distribution, S(D,W,T'), of drop diameter for all values of T' 20 
and W20.5 is shown in figure 14 plotted on normal probability paper. 
Similar plots were made for W20.5, T'222, and W20.8, T' 20; and 
since these distributions were very similar to the one for W20.5, 

T' 20, a statistical test, the chi-square test (reference 17) was 
applied to test the significance of the differences. The results of 
this test showed that the differences between these distributions were 
no greater than would be expected as a result of random sampling; hence 
it was concluded that any dependence of 5 upon Wand T! for values 
of W>0.5 was too small to be reliably indicated by the data. Accord- 
ingly, the distribution S(D,W,T') shown in figure 14 was used for all 
values о? W20.5. | 


In order to define the equiprobability surfaces in the region 
where W is small and D is large, the distributions S'(D,T') 
and R'(W,D,T') must be determined. The distribution, S'(D,T'), of 
the maximum value of D тег icing encounter was examined for T'Z0, 
T! 222, and T'2 32; and the application of the chi-square test to these 
distributions indicated that there was no significant detectable depend- 
ence of 8! on temperature. The distribution, S'(D,T'), for ТІ>0 із 
shown in figure 15. 


In order to determine R'(W,D,T'), the distribution of W for large 
values of D, this distribution was plotted for the following cases: 
D>20, m 20; 5230, T! 20; 1244, T! 20; D220, Т! 222; and 230, T! 222. 
А study of these distributions indicated that there was apparently no 
effect of T' on the variation of R' with W. There was, however, 

a significant effect of D оп the variation of R! with V. 
Figure 16 shows the distributions for Т' 20 and D220, 30, and №. 
Figure 17 shows the parameters defining the Е! distributions Рог 
T! 20 plotted as a function of D. It is likely that the failure 
of the data to show for large values of D a decrease of W with 
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increasing T', as was the case in figure 13, is due to the fact that 
insufficient data are available for large values of D to define the 
function R' when Т! is greater than 22, since it is only in this 
range that the function В becomes strongly dependent upon Т'. Іп 
the absence of more complete information, it has been assumed that 

when T'222, the dependence of Е! on T' is similar to the dependence 
of В on Tt, Accordingly, the values of the ЕБ! distribution param- 
eters (the mode Нова and the fifth percentile Wo...) shown in fig- 
ure 17 are altered to obtain values of the distribution parameters which 
are considered applicable when T' >22. Тһе factor used to alter any 
particular value of W presented in figure 17 is established from fig- 
ure 13 by forming the ratio of the value of W at a selected value 

of T' to the value of М applicable to values of T' $22. The values 
ОҒ Мо вз 219 Wo оқ Which define the function К! were obtained in 
this manner, and are presented in table VII. 


Construction of the equiprobability contours.- After the distribu- 
tion functions have been established, the equiprobability contours Рог. 
various temperature depressions can be determined. The following 
examples illustrate the use of the probabilities Q, R, S, R', and S' 
in determining the equiprobability surface for Pacific coast cumulus 
clouds corresponding to an exceedance probability of 0.01. For a value 
of Pe = 0.01, the equation 


QRS = 0.01 (42) 
applies to large values of W. Similarly, the equation 
QR'S' = 0.01 (A3) 


applies to small values of М. Іп evaluating equations (А2) and (АЗ), 
two equiprobability surface contours are selected. These contours are 
for Т! = 0 and T! = 32. 


Large values of W: Table VIII shows the calculation of W, 
using equation (АР) for various arbitrarily chosen values of D. Com- 
putations are shown when T! = О and when Т! = 32. 


When T' = О, it should be noted that 0 = 1.0. Thus, equa- 
tion (A2) becomes 


В = 0.01/8 (АД) 


The values of 8 necessary in equation (AN) to calculate values of В 
for various values of D are obtained from figure 1%. Values of В 
calculated from equation (Ай) are used to define the desired values 

of W on the curve depicting the distribution function R(W,T') for 

a value of Т! = 0. This function is defined as a straight line on 
Gumbel distribution paper by establishing the values of М for the 
mode and at the fifth percentile. The values, Wo „ga and Wo on used 
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in defining the straight-line distribution of R(W,T') are obtained 
from figure 13 corresponding to a value of Т! = 0. For this partic- 
ular case, the result would be the same as presented in figure 12. 

The desired values of W chosen from the straight-line distribution 
for the calculated values of R are presented in table VIII, part (a). 


When Т! = 32, 0 is found from figure 11 to have the value 
of 0.176. Hence, the equation expressing values of R is 


0.01 
n 0.1768 (45) 
Values of 5 required to solve equation (A5) are obtained again from 
figure 14 for various arbitrarily selected drop diameters. То find 
values of W to which the calculated values of R correspond, it is 
again necessary to establish R(W,T') as a straight-line function 
of W on Gumbel distribution paper. The two points required to 
establish the straight-line variation are the values of W at the 
mode (Мо ва) and at the fifth percentile (Но.ов): These values, for 
а value of Т" = 32, are found from figure 13 to be Wo да = 0.28 
and Wo .o5 = 1.26. The values of W found from the straight-line 
variation of W with R(D,T') are tabulated in table VIII, part (b), 
along with the corresponding values of D, 5, and Е. 


Small values of М: Sections of the equiprobability contours 
corresponding to small values of W are calculated in identically the 
game manner as for large values of W. In both the case of Т! = O 
and of Т! = 32, values of 8! were obtained from figure 15 for cal- 
culating values of R! from the equation 


В! = Р/98! 


The values of W at the mode and the fifth percentile for defining W 
as a function of R'(W,D,T') were obtained from table VII. ‘The calcu- 
lations for Tt = О are presented in table IX, part (а), and the cal- 
culations for Т! = 32 are presented in table IX, part (b). 


Intermediate values of М: From tables VIII and ТХ, equiprobabil- 
ity curves can be plotted for both high and low values of W, as is 
shown in figure 18. For the intermediate values of W, however, the 
final equiprobability contours must be obtained by joining the upper 
and the lower portions by а smooth connecting line, shown by a broken 
line in figure 18. The upper portions of the final contours represent 
maximum values of liquid-water content per icing encounter and the 
corresponding drop diameter; the lower portions represent maximum drop 
size per encounter and the corresponding liquid-water content; and the 
intermediate portions may represent either one of these combinations. 
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The equiprobability contours presented in figure 18 are an 
incomplete set of those presented in figure h(b). 


Саве 5, Layer Clouds in Eastern United States 


The distribution functions.- The distribution of the temperature of 
the most severe icing observation in each icing encounter for this case 
is sHown in figure 19. Since а satisfactory straight-line representa- 
tion could not be found for this frequency distribution, the smooth 
curve on normal probability paper given in figure 19 was used to define 
the function Q(T'). 


The distribution of maximum liquid-water content per icing 
encounter, R(W,T'), was plotted for values of Т' equal to or greater 
than 0, 7, 12, 17, 22, and 27. The curves for T' 20 and T'217 are 
presented in figure 20. The values of Нова and Wo obtained from 
each of the six distribution curves were plotted as “sorctions of Т! 
(fig. 21), and smooth curves were drawn to т" the distribution 
parameters which define the function R(W,T'). 


The distribution of drop diameter for values of W>0.€ and т!>0 
is shown in figure 22. As in case 2, this distribution was found to be 
approximately independent of Wand Т! for values of W greater 
than 0.2. The curve in figure 22 was, therefore, used to represent the 
function S(D,W,T'). 


The distribution of the largest value of D, S'(D,T) observed in 
each icing encounter is shown in figure 23. This distribution curve, 
which was found to be approximately калерендепг of Tt, was used to 
define the function S'(D,T'). 


In order to determine R'(W,D,T'), the distribution of М corre- 
sponding to large values of D, the distribution of W was plotted for 
all values of T! and for values of D equal to or greater 
than 15, 17, 19, 21, 23, and 25. Values of Мо са and Мо оқ for each 
of these distributions were plotted against D as shown in figure 24. 
As a result of the small amount of data available for large values 
of D, there is a considerable scatter in these distribution parameters. 
The smooth curves shown in figure 2) were drawn to represent approxi- 
mately the dependence of R' upon D. To determine the dependence 
of R' upon T', cases of D217 were chosen since this was the 
largest value of D for which sufficient data were available to provide 
a reasonably reliable sample. Distributions of W for D217 and T'20, 
12, and 17 were plotted and the resulting values of Wo.es and Wo „o5 
from these distributions were used with the data from figure 24 to pro- 
vide the basis for the construction of & set of curves defining the 
parameters of the distribution, R'(W,D,T'), (fig. 25). Curves repre- 
senting values of Wo „ag and Wo ов for D217 were drawn through the 
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data points (fig. 25), using the shape of the curves of figure 21 (the 
distribution parameters for Н) ав a guide. The curves for other values 
of D were then drawn from the points determined by the values 

for T' = 0 (obtained from fig. 24) 1п such а way as to conform to the 
pattern established by the curves for D217. These curves define the 
function R'(W,D,T'). 


Values of the probability terms, Q, R, 5, R', and S', were obtained 
with the aid of figures 19, 20, 22, 25, and 23, respectively, and were 
used to construct the equiprobsbility surfaces for values of Pe ! 
from 0.1 to 0.001, which are presented in figure 7. 


Cases 1, 3, and № 


By following procedures similar to those described in the preceding 
examples, equiprobability surfaces alse could be determined for the 
following cases: (1) layer clouds in the Pacific coast area (fig. 3); 
(2) layer clouds in the plateau area (fig. 5); and (3) cumulus clouds in 
the plateau area (fig. 6). The case of cumulus clouds in the eastern 
United States was omitted because of insufficient data. 
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р APPENDIX В 
DETATLS OF THE PREPARATION OF THE PARTIAL 


PROBABILITY CHARTS 


The same cases as were used in describing the construction of the 
equiprobsbility charts will be used to describe in detail the construc- 
tion of the partial probability charts. The general procedure in con- 
structing the charts is first to evaluate the equation 


AP, = Pap Рур Py, (AT) 


The resulting values of AP; for different combinations of liquid-water 
content and intervals of drop diameter and temperature are then arranged 
into a graphical form for convenient use. 


Evaluation of ДР: 


Case 5, layer clouds in eastern United States.~ Values of Pam for 
this case were obtained from figure 19 which provides information as to 
the probability that the temperature depression will be lower than any 
particular value. То find a value of Pap (the probability that the 
т temperature will lie within the interval AT) from figure 19, it is only ` 

necessary to take the difference between the probabilities indicated for 
the two chosen limits of AT. Values of Pan determined from figure 19 
are presented in table X, part (a), for selected temperature intervals. 


Values of Pap (the probability that a drop diameter will lie within 
the interval AD) were obtained in a manner similar to that used in 
obtaining values of Pape Іп this case, however, a composite агор- 
diameter distribution curve was used in order to allow approximately for 
the occurrence of icing encounters in which the conditions of maximum 
severity are determined primarily by drop size rather than liquid-water 
content. The composite distribution curve is presented in figure 26 
where the probability values in this special саве are called 5!!. Ав 
can be seen from figure 26, the composite distribution curve follows 
the distribution of all observed values of drop diameter for the lowest 
90 percent of the distribution &nd the distribution of maximum obseryed 
yalue per encounter (fig. 23) for the highest 1 percent of the cases. 
The intermediate 9 percent of the cases is represented by a straight 
line connecting the two distributions. Values of Pan determined 
through use of figure 26 are presented in table X, part (b), for 
selected drop-diameter intervals. 


4 Іп order to determine values of Fy, (AT), first the distributions 
of maximum liquid-water content per encounter were plotted on Gumbel 
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distribution graph paper for each of the chosen temperature intervals 
‘above -9.59 F in a manner similar to that done in figure 20 for two 
ranges Of values of Т!. Values ОҒ Wo,g4 809 Wo „os, thus could be 
obtained; and these values were plotted ав a function of temperature 
and extrapolated to obtain values of Мо,вз арӣ Wo,os applicable to 
temperature intervals below -9.59 Е. The resulting values of Мо вз 
end Wo.os for all temperature intervals are presented in table X, 
part (с). In addition, there are presented values of W which were 
chosen arbitrarily from the straight-line distribution established by 
the two points, Мова and Wo.o5- Corresponding to the selected values 
of М exist values of Py, (AT) which also сап be read from the straight- 
line distribution defined by the points Wo „ag and Wo „gog These values 
of Py, (AT) are also presented in table X, part (с). 


After the values of Рур» Pap, and Py,(AT) are determined, values 
of АР. can be determined simply by multiplying the constituent prob- 
abilities together for various combinations of liquid-water content and 
intervals of temperature and drop diameter. Such computations are sum- 
marized in table XI for two temperature intervals (32.0° to 20.5° F 
and 20.59 to 10.50 F), two drop-size intervals (0 to 9.5 microns 
and 9.5 to 12.5 microns) and a range in liquid-water-content values. 
The values of АР) listed in the table are ones which сад be obtained 
through use of the partial probability charts. 


Other савев.- The quantities used in calculating values of AP 
for the remaining cases studied were obtained in the same manner as 
ease 5, with one exception. In case 2, the distribution of the greatest 
drop diameter per encounter was used for the upper 5 percent of the 
diameter distribution instead of the upper 1 percent. The reason for 
this deviation was because of the greater frequency of encounters in 
which the most severe icing condition was determined by the maximum 
value of drop size. 


Example Construction of the Partial Probability Charts 


Rather than compute values of AP; every time they are required, 
use of the data presented in tables Х and XI in chart form 18 a great 
convenience. Tn general, the technique of constructing a suitable chart 
for evaluating values of AP; utilizes the relation 


log АР; = log Рур + log Pap + log Py, (AT) (81) 


in such & manner that the chart performs the function of adding the 
right-hand terms in the equation. 


The various steps in the construction of a chart to solve eque- 
tion (Bl) (derived from equation (4)) are shown in figure 27. Fig- 
ure 27(a) establishes the liquid-water content as a function of the sum 
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of the logarithms of Pap and Ру; (АТ). Since the value of Pap corre- 
sponding to a given temperature interval is a constant, the liquid-water- 
content distribution curve is displaced horizontally from the origin by 
the amount of log Pan. The values of Pap, Py, (AT) and Wj used in 
preparing figure 27(a) were obtained from table X, parts (a) and (c), 

for a temperature interval of 32° to 20.5° F. 


The next step in the construction is shown in figure 27(b) where 
the liquid-water-content scale on the right is transformed into а согге- 
sponding scale at the top by use of a line drawn at an angle of 45° on 
the grid system. This line may be labeled the 10-mile horizontal-extent 
line because the ratio of any value on the top scale to the correspond- 
ing value on the right-side scale is unity. For example, in the figure 
а value of Wy = 0.5 on the top scale leads to a value of Wy = 0.5 
on the vertical scale when the lO-mile horizontal-extent line is used. 
For other horizontal extents, however, а factor must be applied to the 
liquid-water content for lO-mile horizontal extent to obtain the value 
appropriate to the extent under consideration. Accordingly, another 
line (derived from reference 3) has been drawn on figure 27(b) to incor- 
porate automatically the factor which should be applied to a 50-mile 
horizontal extent at different values of Pap and Py, (AT). Thus, a 
value of Мі = 0.5 on the top scale would be transformed into а value 
of Wj = 0.82 оп the right-hand scale for a horizontal extent of 
50 miles. The point on the liquid-water distribution curve to which the 
value of W4 = 0.82 applies is labeled point A in figure 27(b). 


The only other factor which must be included in the partial prob- 
ability charts is the probability factor Pap. А method of incorporat- 
ing the term is by using а scale-transformation reference line (similar 
in operation to the horizontal-extent lines). Figure 27(с), the top 
part of which is the same as figure 27(b), is presented to show the 
function of the scale-transformation reference line in determining the 
total sum of log Pam, log Py, (ДТ), and log Pan. Ав an example, a value 
of Wy = 0.5 is taken for which it is desired to find the sum of the 
logs of the three probabilities, Pam, Ру; (AT), and Pap applicable to 
а 10-ші1е horizontal extent, a temperature interval of 32° to 20.59 Е 
and а drop-diameter interval of 0 to 9.5 microns. Accordingly, а value 
of Wi = 0.5 is selected on the liquid-water-content scale and is fol- 
lowed downward, vertically, until an intersection with the 10-mile 
horizontal-extent line is obtained (point 1). From point 1, a line is 
followed horizontally until an intersection with the lLiquid-water-content 
distribution curve for the temperature interval of 329 to 20,59 F 
(point 2). The horizontal distance that point 2 is away from the log 
reference line is a measure of the sum of log Pap and log Py, (AT). If 
а vertical line from point 2 is followed downward until an intersection 
with the scale transformation reference line is reached (point 3), the 
absolute value of the sum of Тов Pam and log Py; (АТ) can be read on the 
scales either to the left or to the right of the point. However, if the 
vertical line followed from point 2 is continued downward until the 
intersection with the drop-diameter interval line is reached (point 4), 
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the sum of the logs of Pam, Ри. (АТ) and Pap can be read from the left- 
hand scale for log of АРІ. The reason why the total can be read is 
that the drop~diameter interval line is displaced vertically downward 
from the scale transformation reference line by the amount of log PAD 
for the diameter interval of 0 to 9.5 microns. The actual value 

of Pap used in the construction was obtained from table X, part (b). 


The final value of log AP; read from the partial probability 
chart in figure 27(c) is the result sought, namely, a solution of equa- 
tion (В1). For most rapid use of the chart, however, it is desirable 
to have values of ДР: obtainable directly. This added convenience 
can be incorporated in a chart by making each integer of the logarithmic 
scales in figure 27(c) correspond to а cycle on semilogarithmic graph 
paper. The partial probability charts presented in figure 10 have been 
constructed to include this feature, and hence will yield values 
of АРА directly. 
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TABLE I.~ CLASSIFICATION OF ICING CONDITIONS 
AND AMOUNT OF DATA AVATLABLE FOR EACH CASE 











бане Geographical 


1 Pacific coast 
2 Pacific coast 
3 Plateau 

n Plateau 

5 Eastern U.S. 
6 Eastern U.S. 











Number of icing 
encounters 


Number of observations 
of liquid-water content 
and drop diameter 


1L 
227 
119 
91 
hol. 
26 


39 
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TABLE ТТ.- INTERVALS OF TEMPERATURE AND DROP DIAMETER 
USED IN THE DETERMINATION OF MARGINAL PROBABILITIES, 
AND MEDIAN VALUES IN EACH INTERVAL TO BE USED IN 

CALCULATING THE MARGINAL LIQUID-WATER CONTENT 


Intervels of temperature, AT Intervals of drop diameter, AD 
(OF) (microns) 


Lower Median of Median of 
limit interval, Ty interval, Di 
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TABLE ІТІ.- COMPARISON OF VALUES OF LIQUID-WATER CONTENT 
FOR PACIFIC COAST CUMULUS CLOUDS, Pe = 0.001, WITH 
‚ CORRESPONDING VALUES FOR INTERMITTENT MAXTMUM 
CONDITIONS LISTED IN NACA TN 1855 












Temperature, T 


(Op) Drop diameter, D 
F 


Liquid or sj ontent, W 
(microns) 


g/m 


TN 1855|Probability analysis 





» 
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TABLE IV.- COMPARISON OF WEIGHTED-AVERAGE VALUES OF 

LIQUID WATER CONTENT FOR LAYER CLOUDS, Pe = 0.001, 

WITH CORRESPONDING VALUES FOR CONTINUOUS MAXIMUM 
CONDITIONS LISTED IN NACA TN 1855 


Liguid—water content, W 


Temperature, T Drop diameter, D (g/m) 


(CF) (microns) 
ron [TN 1855|Probability analysis 
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TABLE V.- CORRESPONDENCE BETWEEN THE EQUIPROBABILITY 
ANALYSIS AND THE ANALYSIS IN NACA TN 1855 














Approximate corre- 
sponding class icing 
condition as defined 
by TN 1855 





















Equipment to 
which condition 
is applicable 


Type of cloud 
formation 






Any part of 
imm|@irplane, such ав 








Tropical cumulus Instantaneous max 


























(т-м) guide vanes, inlet 
Pacific coast ducts, where а 
cumulus None sudden large mass 
Do = = - ~ None of supercooled 
Do ---- Instantaneous normal {water would be 
(I-N) critical, even for 


short duration 


Any critical com- 
ponent of the 






























Pacific coast 


Intermittent maximum 
cumulus 






(II -M) airplane, where ice 
Do = = = < None accretions, even 
Do ---- Intermittent normal though slight and 

(II-N) of short duration, 


could not be 
tolerated 


All components of 
the airplane; that 
is, every part of 













•001 {Continuous maximum 


Layer clouds (III-M) the airplane should 

occurring anywhere „01. None be examined with 

in the U.S.(Pacific the question in 

coast, plateau and | .1 Continuous normal mind, "Will this ‘part 
be affected seri- 


east coast regions) (III-N) 
| ously by accretions 
during continuous 
flight in icing 
conditions 7! 


AExtrapolated from Pacific coast cumulus data. 
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TABLE УТ.- EXAMPLE PARTIAL PROBABILITY CALCULATIONS FOR PACIFIC COAST CUMULUS CLOUDS 


› АТ |32,0 to 20.5120.5 to 10.5 |10.5 to 0.5 10.5 to -9.5 |-9.5 to -19,5|-19.5 to -І0 
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n 
P =) AP; = 0.0012 + 0.0208 + 0.0308 + 0.0116 + 0.0012 + 0.0001 = 0.0657 
1 
NOTE: Values marked with an asterisk, *, denote that the value is actually smaller than the value stated, 


Since such values arise because of termination of the partial probability charts, they are usually 
so small that the values can be neglected in the summation for total probability. 
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TABLE УТТ.- PARAMETERS DEFINING THE FUNCTION, Е!, 
FOR PACIFIC COAST CUMULUS CLOUDS (CASE 2) 


32 to 10 ло ME Ms foe wm 
| т, (FO)| o to 22 12 62 


Liquid. ы content, , 


ore ons 


30 0.70|0.20 0.1810.07 | 0.3210.035| 0.1610.02 
35 .60| .16 | .055 13). O15) ыы 
10 .51| .12 „35 кән eee кете ға 
15 421 .09 „29 тк жек женге жене 
50 .35| .06 „2% кен жетер bea) esee 
55 .28| O45 „19 Жы ЖИЕ 
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TABLE VIII.- EXAMPLE CALCULATIONS OF EQUIPROBABILITY 
CONTOURS FOR SMALL VALUES OF DROP DIAMETER FOR 
PACIFIC COAST CUMULUS CLOUDS (CASE 2) 


(а) Ре = 0.01; Tt = 0; Q = 1.00 


Drop diameter, D [Probability values iLiquid-water content, W 
(microns) (g/m8) 
О 


8 | R | 
1.96 










Drop diameter, D 
(microns) 
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TABLE IX.- EXAMPLE CALCULATIONS OF EQUIPROBABILITY 
CONTOURS FOR LARGE VALUES OF DROP DIAMETER FOR 
PACIFIC COAST CUMULUS CLOUDS (CASE 2) 


(a) Ре = 0.01; T! = 0; Q = 1.00 


— > pes W 


microns 
m a eee T 






„20 . 050 
„12 „083 
307 „143 
204 .250 
„023 „13 
„010 1.00 





(b) Ре = 0.01; T! = 32; Q = 0.176 


у Liguid-water content, W 
Drop diameter, D Probability values (в јаз) 4 
(microns) 


x м 


| 30 0.33 0.172 
32 | „20 „2814 
ho | 


I2 . +7 
47 .057 | 1.00 
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TABLE Х.- VALUES OF THE FACTORS USED IN THE DETERMINATION 
OF PARTIAL PROBABILITIES APPLICABLE TO EASTERN 
LAYER CLOUDS (CASE 5) 


(a) Temperature 


Upper limit of Lower limit of сто s 
Temperature | temperature interval | temperature interval Pat, that 


interval, AT |Value|Probability, Q | Value Probability, Q | temperature 








— сз. — = — s 





(°F) EE of exceeding of exceeding 4... 
AT 

32 to 20.510 1.000 0.480 
20.5 to 10.5 |11,5 520 „345 

| 10.5 to 0.5 |21.5 „175 „112 | 
| 0.5 to -9.5 | 31.5 .063 .017 
| -9.5 to -19.5 |41.5 .016 .013 
j"19.2 to -10.0 [51.5 .003 .003 


(b) Drop diameter 


Drop- Lower limit of Upper limit of 
diameter diameter interval . |diameter interval] | ОСТУ» TAD 


interval, Probability, [Probability,| that drop diameter) 
5'1, of of D, |5'', of lies within the 



















(microns) microndjexceeding D (1! стопа exceeding D | interval AD | 
0 to 9.5 0 0.780 | 0.220 
9.5 to 12.5) 9.5 „480 | • 300 
12.5 to 15.5} 12.5 „230 | „250 
15.5 to 19.51 15.5 .090 | „140 
19.5 to 29.24 19.2 022 | „068 
29.5 to 49.5} 29.5 .0006 | 0214 
19.5 to œ | 49,5 О - 0006 
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TABLE Х.- CONCLUDED 


(c) Liquid-water content 


Liquid -water = 


Temperature content distribu- 


interval, AT 


Liquid—water 


(Ор) tion parameters 
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TABLE XI.- PARTIAL PROBABILITY VALUES FOR EASTERN LAYER 
CLOUDS (CASE 5) FROM DATA PRESENTED IN TABLE X 


“ММ” 













Liquid-water 

content, Wi, 

for AT and AD 
(g/m3) 






Temperature |Diameter 
interval, AT |interval, AD 
(microns) 























«Шортан 2 
/ Y? 22 Т» ә. 
ЕС Lee 






Figure /- Map of the United States Showing approximate boundaries of 
areas used in the geographical classification of icing data. 
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ССС VS АА. ©. < 


Figure 2.— Graphical presentation of the equiprobability surface which 
presents fhe locus of ай combinations of liquid-water content, drop 
diameter and temperature depression below freezing having the 
same probability, е , Of being exceeded in any single icing 

encounter. 
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Е (7'222) 
Е (T's32) 


Liquid-water content, W, g/m? 





Orop diameter 0, microns 


(а) Ра =04 


Figure 3. — Constant-temperature contours defining fhe 
equiprobabllity surface for Case /, (Layer clouds ; 
Pacific coast). 
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(6) Pe г 0.0! 


Figure 3.— Continued. 
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30 0 
Drop diameter, 0, microns 


(c) Pp = 0.00! 


Figure 3.— Concluded. 
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Líquid-water content, И, g/ms 


(а) Ре = 0,1 


Figure 4.— Constant —temperature contours defining the equi- 
probability surface for Cose 2, (Cumulus clouds, Pacific 


coast). 
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Liquid- water content, И, g/m? 
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Figure 4 - Continued. 
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Figure 4. ~ Concluded. 
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А? o F (T =2 
шиау 3212. 
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Liquid -water content, И, g/m? - 


32° F (Т = 0) 
20° Е (T= / 2) 





м 


LLL У И 
О М 


О 
д 8 /6 24 32 


Drop diameter, О, microns 


(а) Ps = 04 


Figure 5. — Constant~temperature contours defining the 
equiprobability surface for Case 3. (Layer clouds, 
Plateau area). 
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(b) P, = 0.0! 


Figure 9. — Continued. | 
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Liquid-water content, W, 97 m3 
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. Figure 5. ¬ Concluded. 
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3220202 F (Т = 0712) 
10°F (T'= 22) 
09Ғ (T'= 32) 
—/0° F (T'=42) 


Liquid -water content, W, g/ms 





0 |. 4/0 20 30 40 
Drop diameter, 0, microns 


(а) P, = O 


Figure 6. — Constant -temperature contours defining the 
equiprob ability surface for Case 4 (Cumulus 
clouds, Plateau area), 
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Figure 6.— Continued. 
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Figure 6. — Concluded. 
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Figure 7. — Constant —'emperature contours defining the 


equiprobability surface for Case 5. (Layer clouds, 
Eastern United States). 
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Figure Z - Continued. 
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Figure 7.— Concluded. 
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Figure 9. — Marginal surface for a particular element of an assumed thermal system. 
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(а) Case 1 — Layer clouds, Pacific Coast. 
Figure 10.— Partial probability charts. 
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(b) Case 2 — Cumulus clouds, Pacific Coast 


Figure 10. — Continued. 
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(с) Case З — Layer clouds, Plateau Region. 


Figure 10.— Continued. 
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| (4) Case 4 — Cumulus clouds, Plateau Region. 
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(6) Case 5 — Layer clouds, Eastern United States. 


Figure 10.— Concluded. 
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Figure 13.— Distribution parameters which define the probability 
function R(W,T) for Pacific Coast cumulus clouds. 
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Figure 14.~ Frequency distribution of drop diameter for observations 
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Figure /'8.—Equiprobability contours constructed 
from the data of tables ZIT and ІХ; 
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Frobability, №, dimensionless 
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Figure 21 — Distribution parameters which define the probability 


85 


function R(W,T') for layer clouds, Eastern United States. 
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Figure 23.— Frequency distribution of the largest value of drop diameter 
observed during each icing encounter; layer clouds, Eastern United 
States. 
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Figure 24.— Distribution parameters defining the function (ШО) 
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Figure. 25.— Distribution parameters defining the function R'(W,D,T'); layer clouds, Eastern 
United States. 
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Figure 26.— Composite frequency distribution of drop diameter used. іп calculating partial probabilities for 
layer clouds in Eastern United States. 
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Liquid-water content distribution 
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(a) Step ! 


Figure 27 — Example construction of а partial 
probability chart. 
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Liquid-water content, W;, grams per cubic meter 
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(b) Step 2 
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Figure 27 — Continued. 
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